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(Dated: August 12, 2008) 

We report on a measurement of the inclusive jet cross section in pp collisions at a center-of-mass 
energy =1.96 TeV using data collected by the DO experiment at the Fermilab Tevatron Collider 
corresponding to an integrated luminosity of 0.70 fb~^. The data cover jet transverse momenta from 
50 GeV to 600 GeV and jet rapidities in the range -2.4 to 2.4. Detailed studies of correlations between 

systematic uncertainties in transverse momentum and rapidity are presented, and the cross section 
measurements are found to be in good agreement with next-to-leading order QCD calculations. 
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PACS numbers: 13.87.Ce,12.38.Qk 



The measurement of the cross section for the inclu- 
sive production of jets in hadron colhsions provides strin- 
gent tests of quantum chromodynamics (QCD). When 
the transverse momentum (pr) of the jet with respect to 
the beam axis is large, contributions from long-distance 
processes are small and the production of jets can be 
calculated in perturbative QCD (pQCD). The inclusive 
jet cross section in pp collisions at large px provides one 
of the most direct probes of physics at small distances. 
In particular, it is directly sensitive to the strong cou- 
pling constant (a^) and the parton distribution functions 
(PDFs) of the proton. Additionally, it can be used to set 
constraints on the internal structure of quarks [H. De- 
viations from pQCD predictions at large pr can indicate 
new physical phenomena not described by the standard 
model of particle physics. A measurement over the widest 
possible rapidity range provides simultaneous sensitivity 
to the PDFs as well as new phenomena expected to pop- 
ulate mainly low rapidities. These data will have a strong 
impact on physics at the CERN Large Hadron Collider 
(LHC) where searches for new particles and higher di- 
mensions suffer from poor knowledge of PDFs [4] . 

In this Letter, we report on a measurement from the 
DO experiment of the inclusive jet cross section in pp 
collisions at a center-of-mass of ^/s = 1.96 TeV. The 
data sample, collected with the DO detector during 2004- 
2005 in Run II of the Fermilab Tevatron Collider, corre- 
sponds to an integrated luminosity of £ = 0.70 fb~^ Q. 
The increased pp center-of-mass energy between Run I 
(^/s = 1.8 TeV) and Run II leads to significant increase 
in the cross section at large pt — a factor of three at 
Pt 550 GeV. The cross section is presented in six bins 
of jet rapidity (y), extending out to |y| = 2.4, as a func- 
tion of jet Pt starting at pt = 50 GeV, and provides 
the largest data set of the inclusive jet spectra at the 
Tevatron with the smallest experimental uncertainties to 
date. The measurement also extends earlier inclusive jet 
cross section measurements by the CDF and DO collab- 
orations 0, Q and improves the systematic uncertainties 
compared to previous measurements by up to a factor of 
two over a range of rapidity up to 2.4 at high px- 

The primary tool for jet detection is the finely seg- 
mented liquid-argon and uranium calorimeter that has 
almost complete solid angular coverage 0]. The cen- 
tral calorimeter (CC) covers the pseudorapidity region 
I77I < 1.1 and the two endcap calorimeters (EC) ex- 
tend the coverage up to \r]\ ^ 4.2. The intercryo- 
stat region (ICR) between the CC and EC contains 
scintillator-based detectors that supplement the coverage 
of the calorimeter. The Run II iterative seed-based cone 
jet algorithm including mid-points 0] with cone radius 
TZ — \J (Ay)2 + (A0)2 = 0.7 in rapidity y and azimuthal 
angle </> is used to cluster energies deposited in calorime- 



ter towers. The same algorithm is used for partons in the 
pQCD calculations. The binning in jet pt is commensu- 
rate with the measured pT resolution. 

Events are required to satisfy jet trigger requirements. 
Only jets above a given pT threshold are kept by the 
highest level trigger (L3). The cross section is corrected 
for jet trigger inefficiencies (always below 2%) determined 
using an independent sample of muon triggered events. 

The jet pt is corrected for the energy response of the 
calorimeter, energy showering in and out the jet cone, 
and additional energy from event pile-up and multiple 
proton interactions. After applying these corrections, the 
jet four momentum is given at the particle level, which 
means that they represent the real energy of the jet made 
out of the stable particles resulting from the hadroniza- 
tion process following the hard pp interaction.. The elec- 
tromagnetic part of the calorimeter is calibrated using 
Z — !• e+e~ events [8]. The jet response for the region 
I77I < 0.4 is determined using the momentum imbalance 
in 7-l~jet events. The pT imbalance in dijet events with 
one jet in jTyj < 0.4 and the other anywhere in 77 is used to 
intercalibrate the jet response in 77, as a function of jet p^- 
Jet energy scale corrections are typically ^--^50% (20%) of 
the jet energy at 50 (400) GeV. Further corrections due 
to the difference in response between quark- and gluon- 
initiated jets are computed using the pythia [9i] event 
generator, passed through a GEANT-based [lH simula- 
tion of the detector response. These corrections amount 
to ~ -j-4% at jet energies of 50 GeV and ~ —2% at 
400 GeV in the CC. The relative uncertainty of the jet pT 
calibration ranges from 1.2% at pT ^ 150 GeV to 1.5% 
at 500 GeV in the CC, and 1.5-2% in the ICR and EC . 

The position of the pp interaction is reconstructed us- 
ing a tracking system consisting of silicon microstrip de- 
tectors and scintillating fibers located inside a solenoidal 
magnetic field of 2 T [6| . The position of the vertex along 
the beamline is required to be within 50 cm of the detec- 
tor center. The signal efficiency of this requirement is 
93.0 ± 0.5%. A requirement is placed on the missing 
transverse energy in the event, computed as the trans- 
verse component of the vector sum of the momenta in 
calorimeter cells, to suppress the cosmic ray background 
and is > 99.5% efficient for signal. Requirements on char- 
acteristics of shower development are used to remove the 
remaining background due to electrons, photons, and de- 
tector noise that mimic jets. The efficiency for these 
requirements is > 99% (> 97.5% in the ICR). After all 
these requirements, the background is < 0.1% in our sam- 
ple. 

The DO detector simulation provides a good descrip- 
tion of jet properties including characteristics of the 
shower development. The correction to the jet cross sec- 
tion for muons and neutrinos, not reconstructed within 



5 



> 

O 

>> 

T3 

H 
Q. 
T3 

~B 

CM 

T3 



10^ 
10® 
10^ 
10^ 
10^ 
10^ 
10 

1 

10-^ 
10-2 
10'^ 
10-^ 
10"^ 
10-® 



D0 Run II 



|y|<0.4 (x32) 
0.4<|y|<0.8 (x16) 
0.8<|y|<1.2 (x8) 
1.2<|y|<1.6 (x4) 
1.6<|y|<2.0 (x2) 
2.0<|y|<2.4 




\Js = 1 .96 TeV 
L = 0.70 fb-^ 

^cone ~ 

— NLOpQCD 

+non-perturbative corrections 

CTEQ6.5M II =ii 



50 60 100 200 300 400 600 

P, (GeV) 



FIG. 1: The inclusive jet cross section as a function of jet 
Pt in six \y\ bins. The data points are multiplied by 2, 4, 
8, 16, and 32 for the bins 1.6 < \y\ < 2.0, 1.2 < \y\ < 1.6, 
0.8 < \y\ < 1.2, 0.4 < \y\ < 0.8, and \y\ < 0.4, respectively. 



jets, is determined using pythia and is 2%, independent 
of Pt and y. The corrections for jet migration between 
bins in pT and y due to finite resolution in energy and 
position are determined in an unfolding procedure, based 
on the experimental and y resolutions. The jet pt res- 
olution is obtained using the px imbalance in dijet events 
and is found to decrease from 13% at px ^ 50 GeV to 
7% at Pt 400 GeV in both the CC and the EC. The 
resolution in the ICR is 16% at pt ~ 50 GeV decreas- 
ing to 11% at Pt ~ 400 GeV. The method to unfold 
the data uses a four-parameter ansatz function [lo| to 
parametrize the pt dependence of the jet cross section 
convoluted with the measured pT resolution and fitted 
to the experimental data. 

The unfolding corrections vary between 20% at a jet 
PT 50 GeV and 40% at 400 GeV in the CC. In the 
EC and the ICR, the corrections are less than 20% at 
Pt ~ 50 GeV, but increase to 80% at the largest pt and 
y. Bin sizes in pT and y are chosen to minimize migration 
corrections due to the experimental resolution. The y 
resolution is better than 0.05 (0.01) for jets with pT ^ 
50 GeV (400 GeV), and leads to a migration correction 
less than 2% in most bins, and 10% in the highest y bin. 

The results of the inclusive jet cross section measure- 
ment corrected to the particle level are displayed in Fig. 1 
in six \y\ bins as a function of pp. The cross section ex- 
tends over more than eight orders of magnitude from 
PT = 50 GeV to PT > 600 GeV. Perturbative QCD 
predictions to next-to-leading order (NLO) in as, com- 
puted using the fastNLO program [l^ (based on NLO- 



JET-t-h [13|) and the PDFs from CTEQ6.5M '14|, are 
compared to the data. The renormalization and factor- 
ization scales {fin and fip) are set to the individual jet 
Pt- The theoretical uncertainty, determined by chang- 
ing fift and fip between pt/2 and 2pT, is of the or- 
der of 10% in all bins. The predictions are corrected 
for non-perturbativc contributions due to the underlying 
event and hadronization computed by pythia with the 
CTEQ6.5M PDFs, the QW tune [Hi, and the two-loop 
formula for as- These non-perturbative corrections to 
theory extend from +10% to +20% at pt ~ 50 GeV be- 
tween \y\ < 0.4 and 2.0 < \y\ < 2.4. The corrections are 
of order +5% for pt ^ 100 GeV, and smaller than +2% 
above 200 GeV. 

The ratio of the data to the theory is shown in Fig. 2. 
The dashed lines show the uncertainties due to the differ- 
ent PDFs coming from the CTEQ6.5 parametrizations. 
The predictions from MRST2004 ^ are displayed by 
the large dashed line. In all y regions, the predictions 
agree well with the data. There is a tendency for the data 
to be lower than the central CTEQ prediction — partic- 
ularly at very large pt — but they lie mostly within the 
CTEQ PDF uncertainty band. The pT dependence of 
the data is well reproduced by the MRST parametriza- 
tion whose systematic uncertainty is slightly smaller than 
that from the CTEQ parametrization. The experimental 
systematic uncertainty is comparable to the PDF uncer- 
tainties. The theoretical scale uncertainty, obtained by 
varying the factorization and renormalization scales be- 
tween — fj,p — pt/2 and /i/? = fip = 2pT, is typically 
10-15%. In most bins, the experimental uncertainties are 
of the same order as the theoretical uncertainties. Tables 
of the cross sections together with their uncertainties are 
given in Ref. 17| . 

Correlations between systematic uncertainties are 
studied in detail to increase the value of these data in 
future PDF fits ^Y?\ and their impact on LHC physics in 
particular. Point-to-point correlations in pT and y are 
provided for the 24 sources of systematic uncertainty. 
The relative uncertainties in percent on the cross sec- 
tion measurement are shown in Fig. 3 for the five most 
significant sources of systematic uncertainty in |?;| < 0.4 
and 2.0 < \y\ < 2.4. The luminosity uncertainty of 6.1%, 
fully correlated in pT and y, is not displayed in Fig. 3. 
The other y bins have similar correlations in shape and 
values between these two extreme bins. The total uncor- 
related uncertainty is < 3% in the CC, and < 15% in the 
EC. 

The two largest systematic uncertainties are due to the 
electromagnetic energy scale obtained from Z — > e'^e~ 
events 0], and the photon energy scale in the CC ob- 
tained using the difference in the calorimeter response 
between photons and electrons in the detector simulation. 
The uncertainty on the photon energy scale is mainly due 
to the limited knowledge of the amount of dead material 
in front of the calorimeter and from the physics modeling 
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FIG. 2: Measured data divided by theory for the inclusive jet cross section as a function of jet pr in six \y\ bins. The data 
systematic uncertainties are displayed by the full shaded band. NLO pQCD calculations, with renormalization and factorization 
scales set to jet pr using the CTEQ6.5M PDFs and including non-perturbative corrections, are compared to the data. The 
CTEQ6.5 PDF uncertainties are shown as small dashed lines and the predictions with MRST2004 PDFs as large dashed lines. 
The theoretical scale uncertainty, obtained by varying the factorization and renormalization scales between = hf = Pt/'2 
and = jiF = '2pT, is typically 10-15% 



of electromagnetic showers in the GEANT-based sim- 
ulation. These two contributions to the jet cross section 
uncertainty are ~ 5% in the CC and 5 — 15% in the EC. 

The large-pT extrapolation of jet energy scale is de- 
termined using the detector simulation with the single- 
pion response tuned to 7-l-jet data. The uncertainty 
rises to 12% (30%) in the CC (EC), and is dominated 
by the uncertainty in the jet fragmentation, estimated 
by comparing the fragmentation models in pythia and 
HERWIG flS^. The uncertainty in 77 intercalibration cor- 
responds to systematic uncertainties associated with the 
procedure to equalize the calorimeter response in differ- 
ent regions of rj in dijet events. These systematic uncer- 
tainties are negligible in the CC because the 77 dependent 
response is calibrated with respect to the CC, but ex- 
tend up to 25% in the EC. Finally, systematic uncertain- 
ties associated with showering effects, due primarily to 
the modeling of the hadronic shower development in the 
detector and differences between pythia and herwig, 
range from 3% at low pT to 7% (15%) at large px in the 
CC (EC). 

To show the potential impact of using point-to-point 
uncertainty correlations in jet pr and y on PDF deter- 
mination, we give in Fig. 3 the uncorrelated and total 
systematic uncertainties as a function of jet pr as a per- 
centage of the jet cross section measurement. The total 
uncorrelated uncertainties are less than 15% and 25% 
of the full uncertainties in the CC and EC respectively. 
The full systematic uncertainties are similar in size to 
the PDF uncertainties (Fig. 2) and the detailed analysis 
of the correlations which have been performed will make 
it possible to further constrain the PDFs. Knowledge of 



these correlations is especially important for constraining 
the PDFs in NNLO pQCD fits where the uncertainties 
due to the dependence on the choice of the renormaliza- 
tion and factorization scales are smaller. The point-to- 
point correlations for the 24 different sources of system- 
atic uncertainties are given in Ref. 



In conclusion, the measured inclusive jet cross section 
corrected for experimental effects to the particle level in 
pp collisions at -^s = 1.96 TeV with C = 0.70 fb"^ is pre- 
sented for six \y\ bins as a function of jet px, substantially 
extending the kinematic reach and improving the preci- 
sion of existing inclusive jet measurements. NLO pQCD 
calculations with CTEQ6.5M or MRST2004 PDFs agree 
with the data and favor the lower edge of the CTEQ6.5 
PDF uncertainty band at large pr and the shape of the 
Pt dependence for MRST2004. A full analysis of corre- 
lations between sources of systematic uncertainty is per- 
formed, increasing the potential impact of these data in 
global PDF fits and on new phenomena searches at the 
LHC. 
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